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Abstract 

We calculate properties of the ground and excited states of nuclei in 
the nobelium region for proton and neutron numbers of 92 < Z < 104 
and 144 < N < 156, respectively. We use three different energy-density- 
functional (EDF) approaches, based on covariant, Skyrme, and Gogny func¬ 
tionals, each with two different parameter sets. A comparative analysis of 
the results obtained for quasiparticle spectra, odd-even and two-particle mass 
staggering, and moments of inertia allows us to identify single-particle and 
shell effects that are characteristic to these different models and to illustrate 
possible systematic uncertainties related to using the EDF modelling. 
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1. Introduction 


Recent experimental studies of nuclei in the nobelium region provided 
rich spectroscopic data 0,3, which, in principle, can be used as a bench¬ 
mark information for extrapolations into the region of superheavy nuclei. 
Numerous theoretical studies are aimed at modelling of these spectroscopic 

data 01,1. @, 0,0,0,0, EH, EH, EH, El, EH, EH, EH, EH, EH, EH, EH so as to 


make such extrapolations as reliable as possible. 

The estimation of theoretical uncertainties is one of the most essential 
aspects of extrapolating nuclear models into exotic nuclei 22], One, fairly 


easy part of it, is the evaluation of statistical uncertainties of observables that 
are related to the uncertainties of model parameters adjusted, in one way or 
another, to experimental data. Another one, very difficult, pertains to those 
systematic uncertainties related to the definition and contents of the different 
terms that make up the models themselves. An obvious strategy, which, 
anyhow, gives us only a limited glimpse on possible systematic uncertainties, 
is to study a set of variants of a given model, and to analyze differences 
obtained for calculated observables. 

In the present study, we aim at such an analysis of results obtained within 
three fairly different energy-density-functional (EDF) approaches. Namely, 
we employ the covariant EDFs 0], with one classic (NL1 EH]) and one recent 
(NL3* EH]) parameter set, Skyrme EDFs 20], with one classic (SLy4 EH]) 
and one recent (UNEDF2 127111 parameter set, as well as Gogny EDFs [281 ]. 
again with one classic (D1S[29|) and one recent (DIM [30;]) parameter set. 

Our goal is thus to determine, present, and compare results obtained 
within these six models for a common set of calculated observables. We 
aim at performing these analyses within the most similar and/or equiva¬ 
lent conditions, so as to meaningfully discuss general qualitative similari¬ 
ties and differences. In all cases, pairing correlations are treated on the 
Bogoliubov level. In the literature, it is customary to label such calcula¬ 
tions as Hartree-Fock-Bogoliubov (HFB) in the non-relativistic cases and as 
relativistic Hartree-Bogoliubov (RHB) for the specific variant of relativistic 
mean-field model used here, thereby emphasizing that HFB-like equations 
are solved instead of simpler HF+BCS equations. Only in the calculations 
with a Gogny force, however, the same effective interaction is used to de¬ 
termine direct, exchange, and pairing matrix elements. By contrast, in case 
of Skyrme EDF and the relativistic approach particle-hole and pairing ma¬ 
trix elements relate to different effective interactions, which can be used to 
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simplify their form and phenomenological adjustment. Also, in the RHB ap¬ 
proach all exchange terms are neglected, whereas in the case of many Skyrme 
parameterizations some are modified in order to improve the description of 
data j25| . None of these formal differences is relevant for our discussion, and 
we will use whenever possible the generic notion of an EDF method for all 
three approaches. Indeed, in all cases the total energy can be cast into the 
form of a functional of normal and anomalous one-body density matrices from 
which the equations-of-motion are then derived by variation. We do not at¬ 
tribute too much of an importance to quantitative similarities and differences 
between the obtained results, especially when the models are compared to 
experimental data. Indeed, a detailed agreement with the data may crucially 
depend on specific model-parameter adjustments, or on various corrections 
taken into account or disregarded. The phenomenological EDF used include 
a limited set of parameters (typically between seven and twenty) and aim at 
a global description of a wide variety of nuclear properties and therefore a 
perfect agreement with experimental data is out of reach at present. Cer¬ 
tainly, in the future all EDF approaches will be improved; here we only look 
into generic properties obtained for selected current global parameterisations 
thereof. 

In the present analysis, we systematically calculated the ground states 
of even-even and odd-mass nuclei from uranium (Z = 92) to rutherfordium 
(Z = 104) and for neutron numbers between N = 144 and 156. The selection 
of this region of heaviest actinides/lightest superheavy nuclei is guided by the 
need for reliable experimental data on spectroscopic properties (in particular, 
on the single-particle energies of deformed one-quasiparticle states) based 
on which the extrapolability of a given theory/functional towards region of 
super heavy nuclei may be judged. In addition, we determined low-lying 
quasiparticle spectra of odd-mass nuclei and low-spin moments of inertia 
of even-even nuclei. The main thrust of the analysis was on the attempt 
to identify single-particle and shell-structure properties of these nuclei by 
looking at many-body observables such as masses, odd-even and two-particle 
mass staggering, and excitation energies. 

The paper is organised as follows. Selected theoretical aspects of our 
calculations are presented in section [2] with four subsections discussing the 
methods related to obtaining results for the Skyrme EDF SLy4 (12.11) . Skyrme 
EDF UNEDF2 (12.21) . Gogny EDFs (I2.3[) . and covariant EDFs (12.4p . The 
results of the calculations are given in section |3j with subsections devoted to 
the Nilsson diagrams (13. ip . quasiparticle spectra (I3.2l and l3.3p . odd-even and 
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two-particle mass staggering (13. 4 1) . and moments of inertia (13. 5|h Conclusions 
are presented in section 31 


2. Calculation details 


2.1. Skyrme energy density functional SLyf 

For the Skyrme EDF SLy4 [26j, the calculations were carried out with 
the Skyrme HFB solver CR8 whose development over the years has been 


documented in Refs. |31|, [32|, |33|, |34j. It uses a 3D coordinate space mesh 


representation of single-particle states along the lines of the solver EV8 


described in Ref. 35j, but is extended in such a way that intrinsic time- 


reversal invariance can be broken and that HFB equations are solved instead 
of HF+BCS. Single-particle states are represented in a cubic box of 32 3 fm 3 
with a step size of 0.8 fm between discretization points. Imposing triaxial 
symmetry, only 1/8 of the box has to be represented numerically, meaning 
that only a 20 x 20 x 20 mesh is to be treated. 

When calculated with SLy4, the ground states of even-even nuclei con¬ 
sidered here are all axial, and the blocked states of odd-A nuclei also re¬ 
main almost axial. All blocked calculations were initialized with the ground 
states of adjacent even-even nuclei. Self-consistent blocking was performed 
by considering the quasiparticle state dominated by a given eigenstate of the 
single-particle Hamiltonian and by exchanging the corresponding columns of 
the HFB U and V matrices after the diagonalization of the HFB Hamilto¬ 
nian, which in turn was constructed using the mean fjelds_of the blocked 
solution from the previous iteration, see, e.g., Refs. [ 


36, 3H 


To avoid mixing of quasiparticle states with different average values of 
the angular momentum component (jy) parallel to the symmetry axis of 
the initial configuration in the diagonalization of the HFB Hamiltonian, the 
many-body expectation value of (Jy) was held fixed with a cranking con¬ 
straint at the value of (Jy) equal to the one of the blocked quasiparticle 
state. As the code CR8 allows for triaxial shapes the mixing cannot be fully 
suppressed. As a consequence, the blocked HFB states are not necessarily 
orthogonal even when they have different average value of (Jy). 

In this respect, blocking the (jy) = 1/2 levels presents a particular dif¬ 
ficulty. Without a cranking constraint, the code CR8 very often converges 
toward a solution where the many-body expectation value (Jy) is close to 
zero, and where the blocked quasiparticle in the spectrum of eigenstates of 
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the HFB Hamiltonian is mixed with other low-lying quasiparticles of differ¬ 
ent (hi). For these states, using or not using the cranking constraint might 
in some cases make a difference of the order of 100 to 200 keV. An example is 
the ground state of 251 Cf. In blocked calculations with cranking constraint, 
there is a low lying (J||) 7r = l/2 + level at 45keV excitation energy above the 
calculated (J||) 7r = 3/2 + ground state. In the calculations without cranking 
constraint, the energy of the 3/2 + state does not change much, but the l/2 + 
level is lowered by about 180 keV and becomes the ground state, in agree¬ 
ment with experiment, but at the expense of the blocked quasiparticle being 
a strong mixture of (j||) = 1/2 and (j||) = 3/2, and of having an angular 
momentum (Jy) that cannot be easily interpreted within the strong coupling 
model anymore. For blocked states with higher (j||), the mixing is always 
much smaller. 

For the coupling constants of the so-called time-odd terms that contribute 
to cranked and blocked states, the same ” hybrid” choice was made as in 

many earlier calculations @, 0, 0 0, 0, 0, 0, 0, 0, 0, 0. They 

were set to the ” native” values dictated by a density-dependent Skyrme two- 
body force for all terms except for those that multiply terms that couple 
two derivatives and two Pauli spin matrices. The latter were set to zero for 

for 


reasons of Galilean invariance and internal consistency, cf. Refs. [31], |38 
further discussion. 

In the present study, neutron and proton surface pairing interactions were 
used, with the strengths adjusted to the three-point, gaps A® of 247 Cf (N = 
149) and Ap 3) of 249 Bk (Z = 97), as defined by Eqns. (J7J) and ([6]) below, 
leading to 


V n = —1240 MeV fm 3 , 
V p = -1575 MeV fm 3 . 


( 1 ) 

( 2 ) 


The pairing-active spaces were limited by the soft cutoffs of 5 MeV above 


and below the neutron and proton Fermi energies as described in Ref. 45 


With that, the proton pairing strength is significantly larger than that given 
by the standard value of V n = V v = —1250 MeV fm 3 used for heavy nuclei 


in previous publications |39|, [7|, |40|, |44j, [42 


01 (which all have been carried 
out using the Lipkin-Nogami (LN) scheme, though). These previous values 
were obtained by adjusting moments of inertia in rotational bands of heavy 
nuclei 0]. There is, however, also some published work on very heavy nuclei 
4a, |46] where the pairing interaction of volume type was adjusted to some 
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A g values of odd-A actinide nuclei. 

In principle, the like-particle pairing interaction should be of isovector 
type, which implies V n — V p . Differences between the adjusted neutron and 
proton pairing strengths can have many reasons: (i) the compensation of 
imperfections of the calculated single-particle spectra, (ii) the compensations 
of the cutoff energy that is chosen to be the same in neutron and proton phase 
spaces (iii) the compensation of the imperfections of the chosen form of the 
pairing interaction itself. Note that if there were none of the above mentioned 
deficiencies, the proton pairing strength would have to be smaller than the 
neutron pairing strength, so as to compensate for the absence of Coulomb 
pairing in our calculations. 

Our attempts to adjust the neutron pairing strength to Ain 251 Cf 
(N = 153), as done in Ref. j2l| and in the present work for UNEDF2, led 
to the values of V n — —1060 and V p = —1565MeVfm 3 . With such much 
weaker neutron pairing strength, in ground states of many of the lighter 
odd-A nuclei, the pairing disappeared. This seems to be connected to an 
anomaly of the calculated single-particle spectra of heavy A = 153 isotones 
that translates into a much larger values of A® for these nuclei than for the 
neighbouring ones. 


2.2. Skyrme energy density functional UNEDF2 

For the Skyrme EDF UNEDF2 j27|, the calculations were performed using 
the symmetry-unrestricted code hfodd (v268h) J47|. The HFB equations 
were solved by expanding single-particle wave functions on 680 deformed 
harmonic-oscillator (HO) basis states, with HO frequencies of Ml x = Ml y = 
8.4826549 MeV and MR = 6.4653456 MeV. This corresponds to including 
the HO basis states up to N x = N y = 13 and N z = 14. 

The time-odd coupling constants of the Skyrme EDF UNEDF2 were de¬ 


termined by the local-gauge-invariance arguments, as defined in Ref. 48 


The HFB equation adopted a quasiparticle cutoff energy of 60 MeV. Within 
this model space, the neutron and proton pairing strengths were readjusted 


/o\ 

to match the experimental values of A" 
tively. This resulted in values of 


n and A p 


( 3 ) 


m 


251 


Cf and 249 Bk, respec- 


V n = -233.889 MeV fm 3 , (3) 

V p = -280.330 MeV fm 3 (4) 

that are larger by about 25 and 50 MeV fm 3 , respectively, than those corre¬ 
sponding to the original UNEDF2 values, where the LN pairing corrections 
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were included, though. 


2.3. Gogny energy density functionals 

Odd mass nuclei were described using the HFB method with full blocking 
and taking into account all possible time-odd fields coming from the poten¬ 
tial part of the Gogny interaction. The original formulation of the density- 
dependent part of the interaction is preserved and therefore no additional 
time-odd density-dependent terms are added to it. This assumption has to 
be verified in the future as there are indications that such kind of terms coiild 


be required to recover basic properties of odd-odd systems [49 


The computer code ATB was used in the calculations 50]. Axial symmetry 


was preserved in the calculations but reflection symmetry was not. Therefore, 
the projection of the angular momentum along the intrinsic axial-symmetry 
axis, = (j’ 11 ), is a good quantum number, but this is not necessarily the case 
for parity. However, many of the states analyzed show a mean value of the 
parity operator, which is rather close to either plus or minus one, allowing 
the labelling of those states with a definite parity. As is a good quantum 
number, states with different values of are automatically orthogonal. States 
with the same value, same parity and similar deformation parameters are 
imposed to be orthogonal to each other by using the traditional technique 
of constraints. The way the constraints are handled does not depend on the 
nature of the constrained operator as only its mean value and gradient are 
required [51] . The code expands the quasiparticle operators in a harmonic 
oscillator basis with 15 shells (680x2 states) although axial symmetry reduces 
the maximum size of matrices to 120 x 2, corresponding to the 0 = 1/2 block. 

In the odd-mass (odd number parity) case, the HFB iterative minimiza¬ 
tion process requires a starting mean-field state with given characteristics. 
An HFB mean-held state with even number parity (obtained with a time-even 
code) was used to generate this starting configuration. The mean values of 
neutron N and proton Z number operators were constrained to the required 
values for the odd nucleus under consideration. To avoid spurious good-parity 
solutions, which could result from the propagation of self-consistent symme¬ 
tries in the HFB method, a small octupole moment was also constrained. 

Ten quasiparticle states with the lowest one-quasiparticle energies coming 
out of this calculation were considered as starting configurations for the time- 
odd code atb. They were obtained with the standard “blocking procedure” 
consisting of swapping the appropriate columns in the U and V matrices 
of the Bogoliubov transformation. Each of the starting configurations was 






labelled with the corresponding G quantum number, which was preserved all 
along the minimization process. The first starting configuration with a given 
value of usually converges to the lowest-energy band head with angular 
momentum I = Q. For the same value of f2, the second starting configu¬ 
ration usually converges to the same first solution, unless an orthogonality 
constraint to that state is imposed (including both the and — states). 
For the third starting configuration for a given value, the orthogonality 
constraints with respect to the first and second states are required. In prin¬ 
ciple, for the rt-th starting configuration 2 {n — 1) orthogonality constraints 
to the states previously obtained are required. Although the number of con¬ 
straints can be relatively high, this is not a formal problem for the gradient 
method used in the code, as handling of constraints is very simple to im¬ 
plement j5l]]. The orthogonality constraint is also important to ensure the 
orthogonality of the physical states. If such orthogonality is not imposed, 
even if the states are different, the relative excitation energies can be shifted. 
In the present calculations only two starting configurations for each value of 
f] were considered. 

For the calculation of the properties of even-even nuclei required in the 
evaluation of the pairing gap parameters, the code HFBaxial was used. It 
has the capability to break reflection symmetry although the phenomenon 
was not relevant in the present calculation. Moments of inertia were deter¬ 
mined using cranked HFB wave functions obtained with the code HFBtri. 
In this code, the axial and time-reversal symmetries were allowed to be bro¬ 
ken in the minimization process. 


2-4- Covariant energy density functionals 

The RHB equations fhi, 8j were solved in the basis of an anisotropic three- 
dimensional harmonic oscillator in Cartesian coordinates. For all nuclei and 
states determined in this work, the same basis deformation of /?o = 0.3, 
7 = 0° and oscillator frequency of huo = 41A -1 / 3 MeV have been used. All 
fermionic and bosonic states belonging to the shells up to Np = 14 and Nb = 
20 were taken into account when performing diagonalization of the Dirac 
equation and matrix inversion of the Klein-Gordon equations, respectively. 
As follows from detailed analysis of Refs. , this truncation of basis 

provides sufficient accuracy of the calculations. 

As the effective interaction in the particle-particle (pp) channel, the cen¬ 
tral part of the non-relativistic Gogny finite-range interaction was used. The 
clear advantage of such a pairing force is that it provides an automatic cutoff 
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of high-momentum components. The motivation for such an approach to the 
description of pairing was given in Ref. j52j. The D1S parametrization of 
the Gogny force was used here. No specific adjustment of its strength was 
used, because it provided a reasonable description of the pairing indicators 
in 249 Bk and 249 > 251 Cf and moments of inertia in 252 ’ 254 No js]. 

Two different covariant EDFs, namely, NL1 [ 23 ] - fitted to the nuclei in 
the valley of beta-stability, and NL3* j24| - tailored towards the description of 
neutron-rich nuclei, were used in the current study. The covariant EDF NL1 
was extensively used in the calculations of rotational bands across the nuclear 
chart (see Ref. j2o[). Covariant EDF NL3* was less tested than NL1 with 
respect to the description of rotating nuclei. However, its global performance 


is well established 53]. Note that so far only these two covariant EDFs were 


systematically confronted with experimental data on single-particle states. 
For example, the study of predominantly single-particle states in odd-mass 
nuclei neighbouring to the doubly magic spherical nuclei was performed in 
Ref. [sil ] within the relativistic particle-vibration coupling model employing 
covariant EDF NL3*. In Ref. [13], the first systematic study of the single¬ 
particle spectra in deformed nuclei in rare-earth region and actinides was 
performed with the covariant EDFs NL1 and NL3*. 

ft is interesting that the overall accuracy of the description of the energies 
of deformed one-quasiparticle states [l_3] is slightly better in the old covariant 
EDF NL1 than in the recent functional NL3*. This suggests that the inclu¬ 
sion of extra information on neutron rich nuclei into the fit of the functional 
NL3* may lead to some degradation of the description of single-particle states 
along the valley of beta-stability. Note that these two functionals well repro¬ 
duce deformation properties of ground states of even-even actinides jgl, [ 20 ] 
and indicate that they are axially symmetric. 

A proper description of odd or rotating nuclei implies breaking of the 
time-reversal symmetry of the mean field, which is induced by the unpaired 
nucleon 55] or rotation j56j. As a consequence, time-odd mean fields and 
nucleonic currents, which cause the nuclear magnetism 57] have to be taken 
into account. In the covariant EDF, time-odd mean fields are defined through 
the Lorentz invariance, and thus they do not require additional coupling 
constants. 

The effects of blocking due to the odd particle were included in a fully 
self-consistent way. This was done within the code CRHB, according to 
Refs. fj58j, H [ho[. The blocked orbital was specified by different tags such as 
(i) dominant main oscillator quantum number N of the wave function, (ii) 
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dominant 12 quantum number of the wave function, (iii) particle or hole na¬ 
ture of the blocked orbital, and (iv) position of the state within the specific 
parity /signature/dominant- 7V/dominant-f2. For a given odd-mass nucleus, 
possible blocked configurations were defined from the analysis of calculated 
quasiparticle spectra in neighboring even-even nuclei and the occupation 
probabilities of the single-particle orbitals of interest in these nuclei. 

Note that in the cases when the calculations of odd-mass nuclei were 
performed only for the definition of the A^ indicators (see Sec. 13.41 below), 
we restricted the analysis to 5-6 one-quasiparticle configurations with ex¬ 
pected lowest total energies, so as to properly determine the ground state of 
an odd-mass nucleus. The calculations confirmed the conclusion of the sta¬ 
tistical analysis of Ref. [l3j that absolute majority of the one-quasiparticle 
configurations are axially symmetric. However, some degree of triaxiality 
was obtained in the ^3/2[622], z/l/2[501] and 7rl/2[400] configurations. 


3. Results 


In this section, we present the results obtained and we look for signatures 
of shell effects in the systematics of many-body observables. This is in ad¬ 
dition to analyzing the eigenvalues of mean-held Hamiltonians as a function 
of deformation, that is, the Nilsson diagrams, which are not observables, but 
provide a useful illustration of the underlying single-particle structure. 

We note here that the calculation of many-body observables implies a full 
self-consistency reached for every individual state, that is, for ground and 
excited states. For example, quasiparticle spectra, which we discuss below, 
always result from calculating differences of total energies, determined sepa¬ 
rately for different many-body self-consistent solutions. Because of that, each 
blocked HFB state may have a slightly different quadrupole, hexadecapole, 
or hi ghe r deformation (see, for example, the results of statistical analysis in 


Ref. 13]), which then feeds back to the mean and pairing fields. In this way, 


the deformation-polarization effects, exerted by one-quasiparticle states, are 
fully taken into account. 

In odd-mass nuclei, quasiparticle excitations were obtained by blocking 
the relevant levels when performing the HFB calculations. The spectra were 
obtained by comparing the total energies of different configurations. The 
procedure closely followed that of Refs. [8j, 37 . 
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3.1. Nilsson diagrams 

The Nilsson diagrams, shown in Figs. [[HIl have been obtained by diag¬ 
onalizing the self-consistent mean-field Hamiltonians corresponding to the 
states that were constraint to a sequence of values of the axial mass quadru- 
pole moment 6(i. O = (2 z 2 — x 2 — y 2 ). Then, by using a simple phenomeno¬ 


logical formula [61 


h = 


47T 


3 R 2 A V 16vr 


-Q ~ Q x 0.009 b 


-i 


(5) 


values of the average quadrupole moments were translated into values of 
the Bohr parameters /3 2 . The numerical factor of 0.009 b -1 corresponds to 
R = 1.2 fm H 1 / 3 and A = 254. Note that the Nilsson diagrams obtained 
in this way constitute only an illustration of single-particle properties of the 
nuclei in the nobelium region. Indeed, going in either direction in the chart of 
nuclei, the relative distances between spherical and deformed levels change, 
closing some of the gaps visible in Figs. [T] ED while opening others. Examples 
for such evolution of the shell structure of spherical states and the deformed 
minima can be found in Refs. fiil . liiil. [(ioi. 64]. Small, but clearly visible 
changes of the Nilsson diagrams can already be spotted when just adding a 
few neutrons or protons, compare for example the Nilsson diagram obtained 
with SLy4 for 254 No, Fig. [1] with the one for 250 Fm presented in Refs [43j, 
or the Nilsson diagrams obtained with NL1 and NL3* for 254 No, Fig. El with 
those for 244 Cm presented in Ref. 


20 


As we can see, for all considered EDFs, the overall positions and defor¬ 
mation-dependence of single-particle levels is fairly similar. In particular, 
deformed shell gaps, which appear near ground-state deformations of /3 2 = 
0.3, occur at particle numbers of Z — 98, Z = 104 and N = 150 and/or 
N = 152 in the majority of the functionals. The only exception is the 
functional NL3*, which is characterized by additional gaps at Z — 102 and 
N = 148. Another deformed gap at /3 2 = 0.2 is observed for Z = 110. 
Moreover, significant differences in important details are also visible. The 
proton deformed shell gaps appear consistently above the one at Z — 100 
that is tentatively inferred from the experimental data, see discussion in 
Refs. @, £ . 


Interestingly, although different EDFs show similar deformed neutron 
shell closures, we observe dramatic differences in the shell structure at sphe¬ 
ricity between the Skyrme EDF UNEDF2 and the other ones, as shown 
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. 13 / 2 + [ 606 ] 
\' 7 / 2 “ [743 
' h/p [ 633 ]. 


/3 2 deformation 


Figure 1: Proton (left panels) and neutron (right panels) Nilsson diagrams of 254 No 
obtained for the Skyrme EDF SLy4 (upper panels) and UNEDF2 (lower panels). At 
spherical shapes, the orbitals are labelled with spherical quantum numbers. For SLy4, at 
large deformations, the deformed single-particle orbitals are labelled by the expectation 
values (j||) of the projection of the angular momentum on the axial-symmetry axis. For 
UNEDF2, these orbitals are labelled by the Nilsson labels tt[Nn z A] determined using 
code hfodd. Solid and dashed lines are used for the positive and negative parity states, 
respectively. 
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single-proton energies e; [MeV] single-proton energies e ; [MeV] 




11/2“ [725} 
7/2 + [613] 
3/2“ [761] 
1/2” [761] ' 
1/2“ [770] 
3/2 + [622] 
l/2 + [620] ' 

9/2“ [734] . 
7/2 + [624] 
1/2“ [501] 
5/2 + [622] ' 
13/2 + [606] 
3/2“ [501]’ 
5/2“ [503] ' 
7/2“ [743] 
l/2 + [631] 

l/2 + [860] 
9/2 + [615] 
11/2“ [725] 
7/2 + [613] 
3/2“ [761] 
1/2“ [761] 
3/2 + [622] 
l/2 + [620] 
1/2“ [770] 

9/2“ [734] 
7/2 + [624] 
1/2“ [501] 
5/2 [622] 
/3/2“ [501] 

' 13/2^ [606] 
5/2“ [503] 
7/2“ [743] 
1/?' [631] 


0 2 deformation 


0.1 0.2 0.3 

02 deformation 


Figure 2: Same as in Fig. [T| but for the Gogny EDF D1S and DIM. All results and 
dominant Nilsson labels were determined using the axial-symmetry code HFBaxial. 
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Single-proton energies e. [MeV] Single-proton energies e. [MeV] 



0.1 0.2 
deformation 


0.1 0.2 
P 2 deformation 


Figure 3: Same as in Fig. [T| but for the covariant EDFs NL1 and NL3*. All results and 
dominant Nilsson labels were determined using axially symmetric code RHB. 
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in Fig. |H Compared to SLy4, spherical orbital lj' 15 / 2 - ( 1 * 11 / 2 +) is lowered 
(raised) by about 2 MeV, and thus their relative positions are inverted, re¬ 
sulting in the spherical shell gaps at N = 152 and 170, whereas other EDFs 
predict shell gaps at N = 164. 

The strong rearrangement of spherical neutron shells observed for UN- 
EDF2 as compared to all other EDFs is a consequence of its rather large 
Cq J and C( J coupling constants of the spin-current tensor terms of this pa¬ 
rameterization [27| . In the N ~ 152 region, the inversion of the spherical 
level sequence substantially increases the number of hlled spherical shells 
for which the spin-orbit partner is empty, thereby increasing the size of the 
spin-current terms. In fact, such behaviour is often found at mid-shell for 


parameterizations with large attractive tensor terms 65j, [66 


The relativistic NL1 and NL3* functionals have the unique feature that 
they predict a large spherical N = 138 gap of about 3 MeV that is absent 
in all non-relativistic calculations. As the sequence of spherical subshells is 
different, for NL1 this gap is located between the IA1/2+ and 2p 9 / 2 + levels, 
whereas for NL3* it is found between the 1*11/2+ and 1 . 7 15/2— levels. 


3.2. Quasiparticle spectra in 251 Cf and 249 Bk 

In Fig. [4l calculated spectra of low-lying band-heads in 249 Bk and 251 Cf 
are shown along with the experimental data. In these results, the spectra 
were obtained by individually blocking relevant quasiparticle orbitals and 
then plotting differences of total many-body energies of obtained nucleonic 
configurations with respect of the total energy of the ground state, that 
is, they are not at all equivalent to quasiparticle energies understood as 
eigenvalues of the HFB Hamiltonian. The ground states were identified as 
nucleonic configurations with the lowest total energies. 

Nuclear configurations of deformed odd nuclei (one-quasiparticle config¬ 
urations) were labelled by means of the standard asymptotic quantum num¬ 
bers f2[iVn 2 A] (Nilsson quantum numbers) that correspond to the dominant 
component in the wave function of the blocked quasiparticle state. 

We used the convention of plotting positive (negative) values for the ex¬ 
citation energies of quasiparticle configurations that correspond to blocked 
quasiparticle states having norms of the second HFB components smaller 
(larger) than 1/2 before blocking. In this way, the states that are predomi¬ 
nantly of a particle (hole) character appear above (below) zero energy. More¬ 
over, these states are always plotted relatively to the ground-state; thus the 
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Figure 4: Experimental and calculated quasiparticle spectra in 249 Bk and 251 Cf, see text 
for the convention used here. Experimental data are taken from Ref. We label the 
state with the full Nilsson label of the dominant component of the wave function only if 
the squared amplitude of this component exceeds 50%. The exception is the 1/2[7] state 
which is strongly mixed. However, the cumulative squared amplitude of the components 
of the wave function with TV = 7 in the structure of this state exceeds 90%. Thus, we 
label it only by principal quantum number N and Q. 
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ground state is plotted identically at the value of zero energy. This con¬ 
vention facilitates the comparison of fully self-consistent results with the 
Nilsson diagrams. For experimental quasiparticle configurations, we follow 
the a ssig nments of particle/hole character as presented in Figs. 34 and 35 of 
Ref. E 


In these two odd nuclei, prominent intruder configurations correspond 
to the proton 7/2[633] and neutron 11/2[725] orbitals. In Ref. j2l[, these 
two orbitals were used as benchmark states to adjust strengths of the spin- 
orbit interactions. We see that without such an adjustment, none of the 
studied standard EDFs places them at the right position. The ground-state 
proton 7/2 [633] orbital, which experimentally is almost degenerate with the 
3/2 [521] orbital, for Skyrme and Gogny EDFs appears about 500 keV above 
the ground state and for covariant EDFs about 200 keV below the ground 
state, with the calculated ground states corresponding the 3/2 [521] orbital 
(or 7/2[514] for the covariant EDF NL3*). 

The neutron 11 /2[725] orbital, for covariant, Gogny, and SLy4 EDFs, 
appears too high and for UNEDF2 EDF too low above its experimental po¬ 
sition with respect to the ground-state 1/2[620] orbital. On the one hand, 
one can say that on the absolute scale these deficiencies are not large. On 
the other hand, they may point to slightly incorrect positions of spherical 
intruder orbitals, from which one would like to infer the shell structure of as 
yet not-reached superheavy nuclei. This analysis shows that detailed struc¬ 
ture of very heavy deformed nuclei may depend on extremely fine details 
of the present-day theoretical models, which very well may be far beyond 
any reasonable possibility of adjusting them precisely enough to available 
experimental data. 

Similarly as in the analysis presented in Ref. |2lj], as an attempt to im¬ 
prove the agreement with the experimental values, we have considered varia¬ 
tions of the spin-orbit parameter Wls of the Gogny EDF D1S that could in¬ 
fluence relative positions of intruder states. Increasing Wls from its nominal 
value of 130 MeV fm 5 reduces the excitation energy of the 11/2“ state while 
it increases the excitation energy of the 9/2“ state in 251 Cf. These changes 
improve the agreement with experimental data for larger values of Wis- In 
the 249 Bk case, the 7/2 + goes down in excitation energy as Wls increases, 
while the 5/2 + and 1/2“ levels go up. As in the 251 Cf case, the comparison 
with experiment seems to favor larger values of Wls ■ This, however, has to 
be contrasted with the analysis of the shell structure of heavy spherical nuclei 
like 208 Pb, which usually calls for weaker spiu-orbit interaction jjjl 63 . 
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Figure 5: Experimental spectra j67| of the even-IV Z = 99 isotopes (left panel) and even-Z 
N = 151 isotones (right panel). 


However, it is necessary to recognize that the studies restricted to spin- 
orbit potential may have internal limitations that come from the fact that 
possible deficiencies in the description of the energies of the single-particle 
states emerging from the central potential, as for example those inferred in 
Ref. [40], are ignored. The fact that standard Skyrme functionals provide 
better description of the single-particle states in the Z — 115 nuclei than 
the ones with the strength of spin-orbit interaction adjusted to experimental 
data in nobelium region [64j may be related to such limitation. 

To explore the sensitivity of the results to the amount of pairing corre¬ 
lations in the system, for the Gogny EDFs we performed calculations where 
the pairing strengths of protons and neutrons were multiplied by factors f p 
and f n , respectively. The first noticeable fact is that increasing the neutron 
pairing strength does not influence in a significant way the spectrum of 219 Bk 
(odd Z ) as it also happens when increasing the proton pairing strength in 
251 Cf. Increasing f p reduces the excitation energy of all levels except for 
the 5/2 + state that remains more or less constant. The comparison with 
experimental spectra seems to favor larger proton pairing correlations. In 


19 





the 251 Cf case, all the levels except the lowest 3/2 + decrease their excitation 
energy with increasing pairing strength. As in the 2J9 Bk case, increasing the 
pairing correlations in 251 Cf improves the agreement with experimental spec¬ 
tra. The same behavior was observed during the readjustment of the pairing 
strength for SLy4. However, on should keep in mind that pairing strengths 
are predominantly defined by the odd-even mass staggering, see Sect. 13.41 
below. 


3.3. Quasiparticle spectra in Z = 99 (Es) isotopes and N — 151 isotones 
In the left (right) panels of Figs. [5HE1 we show experimental and calcu¬ 
lated spectra of the even-A" Z = 99 isotopes (even -Z N = 151 isotones). The 
results are presented in the convention discussed in Sec. 13.21 with experimen¬ 
tal assignments taken from Refs. [ 63 , 68 ]. We see that the experimentally 


assigned odd-quasiparticle configurations are theoretically always obtained 
at relevant low excitation energies. However, similarly as in the case of 219 Bk 
and 251 Cf presented above, the details of the level spacing and ordering vary 
from one calculation to another and are not very well reproduced. 

Although in this region of nuclei, the experimental information is richest 
in the particular isotopic and isotonic chains studied here, it is still quite 
scarce, and often experimental assignments of configurations are still only 
tentative j67j. Nevertheless, we can already see several conspicuous experi¬ 
mental trends. 

In protons, we see the 7/2(633] groimd states and 7/2 [514] excited states 
at fairly constant excitation energies of about 400 keV. For the Skyrme EDF 
SLy4, this feature is very well reproduced, with a significant drop of this 
excitation energy predicted a.t N = 156. For the Skyrme EDF UNEDF2, 
the 7/2 [514] level is also obtained above the Fermi level, with the excitation 
energy gradually decreasing already at N = 152. Note, however, that for 
the Skyrme EDFs, in 219 Bk the relative positions of these two levels were 
not very precisely reproduced, so the nice agreement obtained in the Z = 99 
isotopes might be fortuitous. For the Gogny EDFs, this pair of the levels 
is obtained at roughly correct excitation energies, whereas for the covariant 
EDFs, these levels are fairly well degenerate. 

In neutrons, the 9/2 [734] ground states are for all EDFs studied here 
well reproduced, apart from the Skyrme EDF UNEDF2, which gives the 
7/2[613] ground state with the hole-character 9/2[734] orbitals appearing 
at excitation energy of about 400-500 keV. Two particle-character excited 
quasiparticle states, 7/2[613] and 1/2[620], show excitation energies which 
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Figure 6: Same as in Fig. [5] but for the spectra calculated for the Skyrme EDF SLy4 
(upper panels) and UNEDF2 (lower panels). Since for SLy4 only the mean values of the 
single-particle angular momenta (j ||) and ( j 2 ) (and not full Nilsson labels) were available, 
in this case the labels were assigned by analogy with the results obtained for the Gogny 
EDF D1S, which provides very similar spectra. 
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Figure 7: Same as in Fig. [5]but for the spectra calculated for the Gogny EDF D1S (upper 
panels) and DIM (lower panels). 
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Figure 8: Same as in Fig. [5] but for the spectra calculated for the covariant EDF NL1 
(ripper panels) and NL3* (lower panels). 
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increase with mass. They are correctly reproduced for the Skyrme SLy4 and 
covariant EDF’s; however, the increase with mass is too fast. Correct trends 
of these two levels are also determined for the Gogny EDFs. Two other 
hole-character experimental quasiparticle states, 7/2 [624] and 5/2 [622], have 
excitation energies weakly increasing and decreasing with mass, respectively. 
For the Skyrme SLy4 and covariant NL1 EDFs, their excitations energies 
both increase quite rapidly with mass, whereas for the Gogny functionals the 
mass dependence is rather correct. In all models studied here, the energy 
splitting of these two levels is too large as compared to data. 

From the systematics of the one-quasiparticle levels one can draw a num¬ 
ber of interesting conclusions about the significance of the Nilsson diagram of 
single-particle levels for the calculated and observed spectra of one-quasipar- 
ticle levels. Already in schematic models there is no quantitative one-to-one 
correspondence, as the presence of pairing correlations modifies the spectrum 
of low-lying states. In nuclear EDF calculations, there are additional self- 
consistency effects from the separate optimization of each one-quasiparticle 
state. 

As already outlined in Sect. 13.11 the various EDFs do not always agree 
on the size of the energy gaps between spherical subshclls in the Nilsson 
diagrams for 254 No plotted in Figs. [T] E] and [3l sometimes even not on the 
sequence of the levels. For example, the relativistic NL1 and NL3* functionals 
predict a large spherical N = 138 gap of about 3 MeV. For the non-relativistic 
D1S, DIM, and SLy4 EDFs this gap is much smaller and even disappears 
for UNEDF2. With the non-relativistic functionals one finds spherical gaps 
at N — 148 and N = 152 instead. The appearance and size of the N = 
138 gap is controlled by the position of the spherical 1 in/ 2 + level below 
and the 2g 9 / 2 + level above, except for NL3* for which the lj' 15 / 2 - is pulled 
below the 2g 9 / 2 +. The neutron 9/2[734], 5/2[622], and 9/2[615] single-particle 
levels that emerge from the spherical IJ 15 / 2 —, 2g g / 2 +, and in/ 2 + subshclls, 
respectively, are found just above and below the deformed N = 150 and 
N = 152 gaps in the Nilsson diagrams for 254 No. Looking at one-quasiparticle 
spectra of 251 Cf, as plotted in Fig. [2 one finds that for SLy4, D1S, and 
DIM the distance between the hole-character 9/2[734] and particlc-character 
9/2 [615] levels is very satisfactorily described within 200 keV, whereas for 
UNEDF2 it is much too large, and for NL1 and NL3* it is much too small. 
On the other hand, in the N = 151 isotones shown in Figs. [5HB1 the distance 
between these two levels is on average too large for SLy4, D1S, and DIM, 
again much too large for UNEDF2, correct for NL1, and again much too 
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small for NL3*. This indicates that the spherical N = 138 gap of NL3* is 
probably too large. Still, it is unlikely that it has to be made as small as the 
one found for Skyrme and Gogny EDFs. 

By contrast, the distance between the 9/2[T34] and 5/2[622] one-quasi- 
particle levels in 251 Cf and N = 151 isotones, which are connected to the 
1.7 15 / 2 - and 2g 9 / 2 + shells, respectively, is overestimated by almost IMeV for 
SLy4, D1S, and DIM, by a few hundred keV for NL1 and UNEDF2, but quite 
well reproduced by NL3*. This indicates first of all that the pronounced 
deformed N = 150 gap visible in the Nilsson diagrams for SLy4, D1S, DIM, 
and NL1 should be much smaller. In fact, the large deformed gap at IV = 150 
that is predicted by a substantial number of nuclear EDFs has already quite 
often been attributed to be one of the main causes for the disagreement 
between calculated and observed spectroscopic properties of nuclei in the 
A fa 250 mass region [2,(8, 21, 40, 43 . It is remarkable that the UNEDF2 and 


NL3* functionals that both describe well the distance between the 9/2[734] 
and 5/2[622] levels give very different shell structure at spherical shape in the 
Nilsson diagrams of Figs. [Hand [3] For UNEDF2, the spherical lj’ 15 / 2 - shell is 
about 500 keV above the spherical 2g 9 / 2 +, whereas for NL3* it is the other way 
round, the latter being a unique feature among the functionals studied here. 
Note that for NL3* the relative position of these levels is quickly moving with 
particle number: for the slightly lighter 244 Cm the 1 715 / 2 - is already above 
2(79/2+ level like for the other functionals studied here, cf. Ref. ^2o| . 

Altogether, these findings indicate that on cannot expect to find a unique 
one-to-one correspondence between the spectra of one-quasiparticle levels 
of deformed nuclei and spherical single-particle levels in a Nilsson diagram. 
Different shell structure in the Nilsson diagrams might lead to similar one- 
quasiparticle spectra and vice versa. Indeed, the rearrangement, polarization, 
pairing and single-particle-mixing effects when constructing self-consistent 
one-quasiparticles state make the connection quite complex and apparently 
also slightly EDF-dependent. In addition, it should be recalled that for nuclei 
in the A fa 250 region the spherical configuration corresponds to a maximum 
of the deformation energy landscape and therefore should not be associated 
with a physical state. 

Similar cases can be found for proton levels. All Nilsson diagrams for 
254 No plotted in Figs. H [2] and [3] exhibit a substantial spherical Z — 92 gap 
that is located between the 2f 7 / 2 - and 1 h 9 / 2 - orbitals in non-relativistic func¬ 
tionals and between the IA 3/2 and 1 ha / 2 orbitals in covariant functionals. It 


has been pointed out in Refs. [69|, [T0| that there is no indication for such 
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a gap, which is also visible in the single-particle spectra of lighter spherical 


nuclei 63], in the available spectroscopic data for the spherical 21 ”Th and 


218 U nuclei. However, the 7/2 [514], 3/2 [521], 1/2 [521], and 7/2 [633] one- 
quasiparticle levels of 219 Bk and the Es isotopes, which originate from the 
spherical 1 h 9 / 2 , 2/ 7 / 2 and 1 * 13/2 subshells surrounding the spherical Z = 92 
gap, are reasonably well described by all EDFs employed here, see Figs. HHEJ 
As a consequence, a large spherical Z = 92 gap in the Nilsson diagram of 
254 No is not in apparent conflict with the available data for deformed nu¬ 
clei in the A & 250 mass region. In particular the non-relativistic Skyrme 
and Gogny functionals describe well the relative position of the 3/2 [512] and 
7/2 [514] levels within a few hundred keV. It is only for the relativistic NL1 
and NL3* functionals that the spacing between the 7/2 [514] and 7/2(633] 
one-quasiparticle levels becomes slightly too small, cf. Figs. [4] and [HI which 
points to a slight overestimation of the spherical Z — 92 gap for these func¬ 
tionals. Still, the necessary shift of the spherical shells that one would expect 
to correct for the disagreement between calculation and data will not even 
reduce the gap to the size found with non-relativistic EDFs. 

It is necessary to recognize that the present investigation represents one 
of first steps in the direction of understanding of the accuracy of the de¬ 
scription of one-quasiparticle states in deformed nuclei, which, however, goes 
beyond previous attempts by directly comparing different classes of the EDF 
approaches for the same set of experimental data. Based on the current set 
of experimental data similar accuracy is achieved for employed EDF’s. So 
far statistical analysis of the accuracy of the description of one-quasiparticle 
states employing full set of available experimental data on proton and neu¬ 
tron states has been performed only in actinides and only in the framework 
of CDFT using NL1 and NL3* EDF’s in Ref. jl3[. Although many of the 
states are described with acceptable accuracy, for some states the deviation 
of calculated energy from experimental one exceeds 1 MeV. 

Less systematic studies have been performed in the non-relativistic EDF 
approaches. Global survey of the ground state configurations in odd-mass 


nuclei employing three Skyrme functionals has been performed in Ref. 71 


In Skyrme EDF, the spectra of few actinides and of odd-proton Ho nuclei 
have been studied in Refs. [7, 37f. The spectra of selected Rb, Y, and Nb 


nuclei have been studied in axial Gogny EDF in Refs. 72, 73 


These investigations reveal the same two sources of uncertainties 13 


Although the same set of single-particle states appear in the vicinity of the 
Fermi level as in experiment, their relative positions and energies are not 
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always correct. This is first source of uncertainties. The second source of 
uncertainty is related to stretched energy scale in model calculations as com¬ 
pared with experiment which is due to low effective mass of the nucleon at 
the Fermi level. 

While the solution of the first problem can be attempted in the EDF 
framework, the current analysis with different classes of EDF approaches 
suggest that it is not likely to remove all existing problems in the descrip¬ 
tion of the single-particle spectra (see also the discussion in Ref. 0 ]). More 
comprehensive solution, which would also address the second source of uncer¬ 
tainty, would require taking into account particle-vibration coupling which 
will explain existing fragmentation of the single-particle states and possibly 
compress the calculated one-quasiparticle spectra bringing them closer to 
experiment. Combined with respective re-parametrization of the function¬ 
als it may lead to the functionals with better spectroscopic quality. It is 
clear that existing functionals are biased towards bulk properties since either 
no (CDFT) or extremely limited information on single-particle properties 
(Skyrme and Gogny EDF) is used in their fitting protocols. Even in this 
situation the calculated spectra are not far away from experiment which can 
be considered as a success and a good starting point for future development. 


3-4- Odd-even and two-particle mass staggering 

In order to analyze the odd-even mass differences, the three-point pairing 
indicators 0 (staggering parameters), 


Af = i (B{Z + 1, N) + B(Z - 1, N) - 2 B(Z, TV)) (6) 

with odd Z and even N, and 

= ±(B{Z,N + 1) + B{Z,N~1)- 2 B(Z, N)j (7) 


with even Z and odd N , where B(Z,N) is the (positive) binding energy of 
the nucleus, have been plotted in Figs. IU1 and [TUI respectively. In our HFB 
and RHB calculations, increased sizes of the shell gaps can, in principle, be 
seen in the odd-even mass staggering parameters (JU]) and ((TJ) , because the 
density of single-particle states has an immediate bearing on the calculated 
intensity of pairing correlations. This is under assumption that the differ¬ 
ences obtained by blocking different orbitals along the isotopic or isotonic 
chains have a lesser impact on the odd-even mass staggering. Irrespective 
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Proton number Z 

Figure 9: Three-point proton odd-even mass staggering, Eq. ©, shown for the odd-Z and 
even-iV nuclei in the nobelium region. Experimental values are based on the AME2012 
atomic mass evaluation j75|]. 


of the detailed reproduction of the experimental values of these parameters, 
calculated results may thus illustrate the shell structure corresponding to the 
EDFs studied here. 

Experimental results, shown in panels (e) of Figs. [9] and [TO], indicate that 
in the nobelium region, values of the staggering parameters are within the 
range of 500-700 keV. On closer inspection, we see several trends in the mass 
dependence of these parameters, which may indicate variations of the shell 
structure due to the influence of the level density on pairing correlations, or 

/o\ 

due to other fine structural effects. In particular, values of A p ; seem to have 
a small dip for the N = 146 isotones at Z — 95, are fairly constant in the 
N = 148 isotones, and gradually decrease with mass in the N = 150-154 
isotones. None of these values indicate a particularly significant shell-gap 

/o\ 

opening near Z = 100. Similarly, small dips in A„ , which show up in the 
Z = 96-98 isotopes at N = 149 and in the Z = 100 isotopes at N = 153, do 
not point to a particularly large shell gap at N = 152. 

This lack of large variations in odd-even mass staggering is at variance 
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Figure 10: Same as in Fig.[9]but for the neutron odd-even mass staggering, Eq. (0, shown 
for the odd-iV and even -Z nuclei. 


with the analysis of two-particle mass staggering given by quantities 
6$ = 2B(Z,N)-B{Z+2,N)-B(Z-2,N) = S 2p (Z, N)-S 2p (Z+2, N) (8) 
and 

= 2B(Z,N)-B(Z,N+2)-B(Z,N-2) = S 2n (Z, N)-S 2n (Z, N+2), (9) 

which were typically used to identify two-nucleon shell gap s in experiment 
and in calculations for spherical [62.1, 76j] and deformed 77, i>j shell closures 
in the predictions of Skyrme and covariant EDFs. As discussed in Ref. j2|, 
experimental values of o 2 n show clear maxima for the Z = 96-102 isotopes 
at = 152 and those of d\ p exhibit maxima for the N = 148-150 isotones 
at Z = 98 and for the N = 152-154 isotones at Z = 100, see also Figs. [Til 
andfl2l 

When looking at the most pronounced features of the calculated odd- 
even mass staggering shown in Figs. |9] and [TOj we see that minima of 
can be seen at Z = 95 (NL1 and NL3* EDFs) and Z — 97 (DIM, SLy4, 
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Figure 11: Two-proton mass staggering, Eq. ([8]). shown for the even-Z and even -N nuclei 
in the nobelium region. Experimental values are based on the AME2012 atomic mass 
evaluation [75]. 


and UNEDF2 EDFs). For the Gogny and Skyrme EDFs, these minima 
disappear at higher neutron numbers and rather monotonic trends are then 
obtained. Similarly, minima of A„ appear at N — 149 (NL1 and NL3* 
EDFs) or N = 151 (DIM, UNEDF2, and SLy4 EDFs); in the latter case, 
in lighter isotopes they tend to shift to N = 149. For the calculated two- 
proton-staggering indicators ()8l) , covariant EDFs, NL1 and NL3*, exhibit 
very strong maxima at Z = 96, at variance with the data, whereas the non- 
relativitic EDFs, DIM, SLy4, and UNEDF2, reproduce experimental maxima 
at Z = 98 in the N = 146-150 isotones but fail to shift these maxima to 
Z = 100 in heavier isotones. This conspicuous experimental fearure thus 
remains unsolved. The calculated two-neutron-staggering indicators ([§]), do 
not reproduce experimental maxima occurring at N = 152. These results 
illustrate the fact that none of the studied EDFs reproduces the experimental 
trends in shell gaps extracted from the two-particle indicators (JH]) and (J9]). We 
note here that the inclusion of the LN method into the calculations renders 
pairing correlations much less sensitive to the shell structure. Therefore, one 
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Figure 12: Same as in Fig. [TT] but for the two-neutron mass staggering, Eq. ([9|). 


/q\ 

then obtains fairly structureless trends of A [ p and An 1 21], although for 


,(3) 


covariant EDFs, one at the same time obtains a significant improvement of 


the overall agreement with experimental values 20 


3.5. Moments of inertia 

The moments of inertia may constitute another independent indicator 
of the shell structure. This is true under the assumption that their val¬ 
ues, similarly as for the odd-even staggering, are dictated by varying pairing 
correlations. Then, larger shell gaps would induce weaker pairing and thus 
larger moments of inertia. 

The kinematical moments of inertia were calculated as 

J ( i) = (10) 

CJ_L 


where co± is the value of the constrained rotational frequency and ( J ±) is 
the expectation value of the component of angular momentum, both in the 
directions perpendicular to the axial-symmetry axis. The dynamic moment 
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of inertia is not considered due to larger experimental uncertainties conse¬ 
quence of higher order derivatives in its definition. Moreover, contrary to 
the kinematic moment of inertia the dynamic one does not depend on spin. 
Thus, the kinematic moment of inertia provides stricter constrain of model 
predictions. 

All calculations of moments of inertia were performed at hui± = 20 keV. 
In Figs. d3] and [U they are compared to experimental values determined as 


J« = 


3 h 2 
E2+ 


( 11 ) 


In nuclei where at the bottom of rotational bands the experimental 2 + levels 
have not been seen, we used values extrapolated [7£| by the Harris formula. 
In this region of nuclei, the low-spin moments of inertia turn out to be very 
weakly dependent on the angular frequency, and, therefore, a specific method 
of extracting them from experiment is not essential. 

The experimental values shown in panels (e) of Figs. [13] and [14] show clear 
maxima of J (v> in function of Z at Z = 100 in the N = 152 isotonic chain, 
and at Z = 98 in the N = 150 isotonic chain, as well as function of N at 
N = 152 in the Z = 100 isotopic chain. These maxima only partly appear 
in those nuclei that show maxima of the two-particle staggering indicators, 
discussed in Sec. 13.41 

In our theoretical calculations, weak maxima of ./-') are obtained at Z = 
98 for N = 146-150 isotonic chains (DHV@, UNEDF2, and SLy4 EDFs) and a 
stronger maxima at Z = 96 for N = 146-152 isotonic chains (NL1 and NL3* 
EDFs), whereas in heavier isotonic chains we see only a gradual increase, 
without indications of increased shell gaps. Similarly, in the isotopic chains 
maxima appear at N = 148 (NL1 EDF) and merely kinks appear at N = 150 
(DIM, UNEDF2, and SLy4 EDFs) 

Comparing moments of inertia, Figs. [[3] and [14] with the Nilsson di¬ 
agrams, Figs. [THHl we see that our calculations with different models and 
forces seem to exhibit rather nice correspondence between the respective 
proton (neutron) single-particle shell gaps and peaks/kinks in the moments 
of inertia along the isotonic (isotopic) chains. 

In Fig. [131 f° r the covariant EDFs NL1 (a) and NL3* (d), the peaks at 
Z — 96 obtained for the N = 144-152 isotonic chains can be associated with 


1 Results for D1S can be found in Ref. [l3 


32 






80 


70 


60 

50 

> 

® 40 

CM 

•K 

Jr" so 

~^~3 

70 


60 
50 
40 

92 94 96 98 100 102 104 92 94 96 98 100 102 104 92 94 96 98 100 102 104 

Proton number Z 

Figure 13: Kinematical moments of inertia of yrast rotational bands of even-even 
nuclei in the nobelium region plotted in the isotonic chains as functions of Z. Theoretical 
values, calculated at the angular frequency of fiio± = 20keV, Eq. m, are compared to 
experimental values extracted from energies of the 2 + states, Eq. (fTTl) . Experimental 
energies are taken from Ref. [67] except for 252 Fm, which is taken from Ref. [7fjj |. 

the shell gap that in Fig. [3] opens up at Z — 96. The peak moves to Z — 104 
for N = 154 and 156 chains, see Fig. [3j For the Skyrme EDFs SLy4 (c) 
and UNEDF2 (f), the peak/kink at Z = 98 may be associated with a shell 
gap at Z = 98 visible in the Nilsson diagram of Fig. [I] For the Gogny EDF 
DIM, the correspondence is not as clear as that visible in other cases. It 
is to be noted that all our calculations predict that for most of the isotonic 
chains, values of J {V) peak at Z — 104. This corresponds to the proton gap 
at Z = 104 that is clearly visible in Figs. [EH3J 

In Fig. [U for the covariant EDF NL3* (d), values of show pro¬ 
nounced maxima at N — 148 for isotopic chains of Z — 96-104. For NL1 
(a), the maximum becomes a kink occurring at N = 150. Neutron numbers 
of N = 148 and 150 correspond nicely to the shell gaps shown in Fig. |3] 
Non-relativistic functionals predict either a peak at N = 150 or a plateau for 
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Figure 14: Same as in Fig. [CJ1 but plotted along the isotopic chains as functions of N. 


N = 150 and 152, particularly for the Z = 100-104 isotopic chains, shown 
in Fig. [L4j 

It is necessary to recognize that for a proper reproduction of experimental 
moments of inertia, the inclusion of the LN method into the calculations 
appears to be more important in covariant 80] than in non-relativistic EDFs. 


The LN method renders the values of calculated moments of inertia much 
closer to the data 20, 81 and at the same time much less sensitive to the 


underlying shell structure, at variance with the data, cf. the discussion in 


Ref. 21 


4. Conclusions 

In this work we have used three different EDF approaches, each for two 
different parameter sets, so as to gain some insight into the degree of sys¬ 
tematic uncertainties that are related to applying these approaches to spec¬ 
troscopic properties of heavy nuclei at the gateway to superheavy region. 
On the one hand, we concluded that the overall coarse description of several 
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spectroscopic properties of nuclei in this region is, in general, correct. On the 
other hand, we identified numerous smaller or larger differences between the 
results obtained within these different EDFs. In particular, none of the stud¬ 
ied global EDF parameterisations precisely describes variations of the shell 
structure seen in experiment. This can be associated with small deficiencies 
of the obtained deformed shell properties. 

On the scale of very precise spectroscopic experimental data, differences 
between various EDF approaches are large. They are of the same order as 
the degree of differences with experiment. This points to still fairly large 
systematic uncertainties inherent to the models currently in use. Moreover, 
none of the studied models could at present be identified as the one that 
systematically performs significantly better than another one. However, it is 
interesting that on the average the accuracy of the description of the excita¬ 
tion energies of deformed one-quasiparticle states in the nobelium region is 
substantially better than the accuracy of the description of the energies of 
dominant single-particle states at spherical shape in odd-mass nuclei neigh¬ 
bouring to doubly magic 208 Pb (even as compared with the calculations which 
include particle-vibration coupling [lil, 821 . 

The obtained results suggest that further work on improving the perfor¬ 
mance of the EDF methods is very much required. First, one can hope 
that within the existing forms of EDFs, one can still find better global 
parametrizations. For Skyrme EDFs, this route has already been explored 
and a negative conclusion was reached [T3], but for covariant and Gogny 
EDFs similar work has not yet been performed. Second, one can hope that 
various beyond-mean-held corrections, not included in the present analysis, 
may have strong impact on the results, and thus modify the current conclu¬ 
sions. At least that happens in the covariant framework, where accounting 
for the (quasi)particle-vibration coupling improves substantially the accu¬ 
racy of the description of predominantly single-particle states in spherical 


medium and heavy nuclei |54|, [15]. In addition, there is still a possibil¬ 


ity of building new functionals, with beyond-mean-held effects incorporated 
from the very beginning. Again, for the Skyrme EDFs, employed together 
with odd-particle polarization effects included, a negative conclusion has re¬ 
cently been reached j82|]. However, so far in the covariant and Gogny EDF 
approaches no such studies have been performed, and such route has to 
be explored. Third, one can also attempt building new classes of EDFs, 
where systematic expansions within different schemes are used. Only very 
recently this route was started to be explored in non-relativistic EDFs, see, 
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Refs. [83|, |8J, [85], 186], [87], [88|, [89|, [90|, [9l|, so its eventual impact on the physics 


discussed here is not yet known. 
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